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The exact solutions of the first order differential equation with delay are derived. The equation 
has been introduced as a model of traffic flow. The solution describes the traveling cluster of jam, 
which is characterized by Jacobi's elliptic function. We also obtain the family of solutions of such 
type of equations. 
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In this paper, we investigate following differential- 
difference equation, 



— x n (t + t) = f(x n+1 (t) - x n (t)), 



(1) 



where r is a real positive constant called "delay" . The 
index n takes integer. The set of equations of this type 
has been introduced for a car-following model of traffic 
flow HI |H . In that case x n is the position of the nth car. 
Equation (|l|) has been popular in many physical phenom- 
ena of relaxation towards an optimal equilibrium state, 
such as relaxation effect in gases, chemical reactions and 
synchronization problem. 

As a model for traffic flow, f(x) = tanh(x) is the rea- 
sonable choice H d}, which was first introduced in the 
model (|], 



d 2 . , 

T — X n (t) 



f(x n+ i(t) - x n (t)) - ^-x n (t) 



^^-^.^ „„ vvy dfI ,rr ( 2 ) 

This model may have some relation to Eq.(Q) for rela- 
tively small tH. Many studies have been made to the 
model of Eq.(g) with /(x) — tanh(x), which has the sta- 
ble traveling cluster solution as traffic jam || . The travel- 
ing cluster solution is characterized as a soliton of mod- 
ified Korteweg-de Vries (MKdV) equation in the vicin- 
ity of the critical point |(|. One interesting question is 
whether Eq.(Q) has such traveling cluster solution or not. 
Our simulation result suggest the existence of such solu- 
tion, which describes a traveling cluster moving backward 
with the velocity v = l/(2r) 0. Recently, we have an 
information that Igarashi, Itoh and Nakanishi have found 
the exact solution of Eq. (|l]) , which characterized by the 
theta function In this paper we also derive a series 
of analytic solutions of traveling cluster in the context 



of our work, and confirm their stability by the numerical 
simulations. 

It is convenient to introduce new variable h„(t) = 
Xn+i{t) — x n {t) and rewrite Eq.llJ), 



dt 



h n (t 



f(h n+1 (t)) - f(h n (t)) 



(3) 



We start at the linear theory assuming that the ampli- 
tude h n (t) is infinitesimal and f(h n ) = /(0) + f'(0)h n . 
Equation (||) becomes 



dt 



h n (t + t) = h n+ i(t) - h n (t) . 



(4) 



Here we set /'(0) = 1 without loss of generality. We 
investigate the solution of the form 



h n {t) — exp(ian + iuot) 
for real a and u). Inserting Eq. (|5|) to Eq. 

sin a/2 1 
a/2 ~ 2r' 



and 



2ujt 



(5) 

we obtain 
(6) 

(7) 



Equation (g) has a solution if 2t is larger than 1, which 
means r = 1/2 is critical. If r is a little bit larger than 
1/2, Eq.(^) has only two solutions a = ±A. Then, we 
obtain the solution of Eq. (Q) as 



h n (t) = exp±iA(n+ 

It 



(8) 



This represents a traveling wave solution with the veloc- 
ity 1/(2t) in the space of index n, which is treated as a 



1 



continuous variable. The wave moves backward against 
the numbering direction, which appears as the traveling 
wave in the real space moving backward in the flow of 
x n . The above analysis is first given by Whitham 0]. 

Now let us investigate the exact traveling wave solution 
of Eq.(||). We treat the index n as a continuous variable, 
and change the notation h n (t) to h(n,t). We introduce 
new variables on the moving flame of traveling wave as 
u = n + vt, where v is the velocity of the traveling wave. 
We search the solution which does not change its form 
on this flame. We define the amplitude of traveling wave 

H(u) = H(n + vt) = h[n,t). (9) 

Eq.(pl) for the amplitude H(u) is expressed as 



v H(u + vt) = f(H(u + 1)) - f{H{u)). (10) 
au 

Replacing u by u — 1/2, we get more symmetric form 
v^-H(u + a) = f(H(u + 1)) - f(H(u - 1)), (11) 



where 



(12) 



We investigate the solution under the condition a = 0. 
This means that the traveling wave of such solutions 
propagates backward with just the same velocity as the 
linear theory, 



1 

2r* 



(13) 



Now we present the definite form of exact solutions 
giving concrete examples of f{x). First we take f(x) — 
tanh(a;), which is a suitable choice for a model of traffic 
flow. Introducing a new amplitude 



G = f(H), 
Eq.(O) with a — is rewritten as 



(14) 



dG(u)ldu 1, „. 1. , . 

^T^=G(u + -)-G(u--). (15) 

We can easily find a solution of Eq. ( |l5| ) in the form 

G(u)=/3sn(au,k), (16) 

where sn is Jacobi's elliptic function with modulus k. The 
parameter a is determined by 



sn(a/2, k) 
a/2 



1 

2? 



and (3 is given by 



(3 = ±k 



4r 



(17) 



(18) 



Eq.(|l7|) has a real solution only if l/(2r) < 1. In the case 
of k = 0, Eq.(|l7|) reduces to the result of linear theory, 
Eq. (^|) . The modulus k is a free parameter of the solution, 
which indicates the existence of many solutions for the 
same traveling velocity v — l/(2r). The relation between 
the modulus and solutions is discussed later. Next we 
take f(x) = tan(x). In this case Eq. (|l5|) is replaced by 



dG{u)/du 



G(u+-)-G(u--). 



(19) 



1 + G(u) 2 ' 2' T 

It is also easy to see that Eq. (|li]) has a solution 

G(u) = ±k^-cn(au,k), (20) 
4r 

The pa- 



(21) 



where cn is another Jacobi's elliptic function 
rameter a is determined by 

sd(a/2,fc) _ _L_ 
_ 2r' 

which also reduces to the result of linear theory, if we 
take k = 0. 

The above two solutions suggest us to expect the exis- 
tence of another solutions using elliptic functions. Actu- 
ally, we have found the family of solutions of Eqs. (|l5| ) 
or (19). The solutions of this family are denoted in the 
form 



(au, k), 



1 

2t' 



(22) 



(23) 



G(u)=±C k — 
4r 

where a is determined by 

A(a/2,k) _ 
a/2 " 

In the above, g and A are appropriate elliptic functions. 
The solution is given by each set of (g,A,Ck) in Tables || 
and|Ti]. 

TABLE I. Solutions for f(x) = tanh(x) 



g 


sn 


ns 


sc 


cs 


ds 


cd 


dc 


nc 


A 


sn 


sn 


sc 


sc 


sd 


sn 


sn 


sd 


c k 


k 


1 


vT - k 2 


1 


1 


k 


1 


Vl — k 2 






TABLE II. Solutions for /( 


x) = tan(a;) 




g 


cn 


nd 


dn 


sd 


A 




sd 


sc 




sc 




sd 


c k 




k 


y/1 


- k 2 




1 




feVl - k 2 



Thus, all Jacobi's elliptic functions are the solutions 
for either Eq. ( |l5| ) or (|l^). This result is well un- 
derstood from the following fact. Jacobi's elliptic func- 
tions g(cm, k) are connected each other by Jacobi's imag- 
inary transformation: a ia, the imaginary transfor- 
mation of modulus: k — > ik and the translation of u: 



2 



U — ► u + K(k)/a, where K(k) is a quarter of the period 
of elliptic functions. These transformations preserve the 
form of Eq. ( |l5| ) or (|l9|), or exchange each other. Each 
pair of two elliptic functions connected by the translation 
gives an equivalent solution. The relation of solutions is 
shown in Fig.|l|. 



,CS y 




FIG. 1. The relation of elliptic solutions for Eq. ( |lE| ) de- 
noted by — sign and Eq. dl9) denoted by + sign. They 
are connected by transformations k — > ik, a — » ia and 
u — > u + K(k)/a, drawn by thin lines, dashed lines and thick 
lines, respectively. They form a double hexagonal structure. 
Each pair of two elliptic functions connected with a thick line 
gives an equivalent solution. The functions in the left and 
right sides divided by the vertical center line are inverse with 
each other. 

The explicit form for the solution of Eq.^) is given by 
H = f^ 1 (G). In this step some solutions for G needs 
appropriate interpretation for the realistic meaning. The 
solution G ~ cs for G — tanh(iJ) is the case, which has 
the region of u where \G\ > 1. On the other hand, the 
solution G ~ sn has a realistic meaning in the whole 
range of it. 

Equation ( |l5| ) or (19) includes the solutions for an- 
other choice of f(x) besides tanh(a:) and tan(ir). Actu- 
ally, f(x) — coth(x) leads to Eq. (|15|) , and f(x) = cot(x) 
leads to Eq.(19) with v — > —v. Jacobi's elliptic func- 
tions also provide the solutions for these systems. This 
result presents some interesting picture. Let us compare 
the systems controlled by f{x) — tanh(x) and coth(x). 
At first sight, the behavior of these two systems seems 
different from each other. We note sn and ns are the solu- 
tions for both systems, which obey Eq.(p"5|). The solution 
H ~ arccoth(ns) in the system of f(x) — coth(x) has es- 
sentially the same form as the solution H ~ arctanh(sn) 
in the system of f(x) = tanh(a;). As the result, these two 
systems may have the same global phenomenon in spite 
of the different type of local interaction. This analysis is 



not simply applicable to the solution nc in the system of 
f(x) = coth(a;), because cn is the solution of the system 
of f(x) = tan(x), which obeys the other Eq. (|T9|) . 

We discuss the meaning of the modulus in our solu- 
tion. The modulus k determines the period of elliptic 
functions, which is related to the number of traveling 
cluster and the boundary condition. For example, k = 1 
gives kink like solution corresponding to the boundary 
condition G{— oo) = — G(oo) , 



G{z) = tanh(a/2) tanh(cra) 
with a satisfying 

tanh(a/2) _ 1 
a/2 ~ 2r" 



(24) 



(25) 



We perform the simulation to check the stability of 
our analytic solution. Fig.^| shows an elliptic solution for 
f(x) = tanh(a;) given by sn 



H(u,k) = arctanhffc — sn(cra, k)), 
4r 



(26) 



with t = 0.501, k = 0.9965 and a = 0.15522 together 
with the result of simulation for Eq. (0) performed in the 
periodic boundary condition with a suitable initial con- 
dition. 
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FIG. 2. Simulation result of h n in the system of 
f(x) = tanh(a^) is shown by diamonds together with the an- 
alytic solution given by sn. The curve of analytic solution is 
shifted by u — > u+constant. 



We are convinced the set of Eqs. dl5| ) and (19) has 
much more solutions constructed with elliptic functions. 
As a fact, the solution found by Igarashi, Itoh and Nakan- 
ishi Q is rewritten in our formulation as 



sn(q(u + a)) + sn(q(M ^ a)) 

G{u) = [3 — r h7, 

sn(cra) 



(27) 
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where a is a free parameter. Actually this satisfies 
Eq. ( |l5| ) . Their solution includes Eq. ([l6]) as a special case 
if we take aa = K(k)/2 ||. The parameter a is deter- 
mined by 



a cn(cra)dn(aa) x+X- 



4r 



where 



sn(aa) 



x± = 1 



x + — X- 



= 0, 



sn 2 (aa) 



sn 2 (a(a ± ^)) 
Then [3 and 7 are given by 

Q 



= ± 



4rsn(aa) ' 



and 



± 



x+_ 

x+ 



(28) 



(29) 



(30) 



(31) 



This solution represents the asymmetry of the widths of 
upper and lower plateaus for traveling cluster in con- 
trast of the solution given by a single sn shown in FigJ^. 
Some combinations of elliptic functions in the numerator 
and denominator of Eq.(p7f) are solutions, those are con- 
structed by the three transformations represented in Fig. 
|l|. This type of solutions has asymmetry in contrast of 
the solutions given by a single elliptic function in Tables || 
and [f| Probably, the set of Eqs. @ and (|l|) has some 
algebraic structure in the space of solutions constructed 
with elliptic functions. 

We remark all these solutions have the common ve- 
locity v = l/(2r). Numerical simulation shows that the 
traveling cluster solutions of Eq. (|l|) preserve their veloc- 
ity as v — 1/(2t) in the deformation of f(x) beyond 
tanh(rr). This fact suggests v is some invariant quantity 
of the structure in the set of the solutions. 

The set of differential-difference equations for G: Eqs. 
( |l5| ) and (|l9| ) offers rich contents of the system with 
traveling cluster solutions which characterized by elliptic 
functions. Eqs. (|l5|) and (|l9|) are related to some soliton 
systems. Our equations can be derived as the traveling 
wave equations for such soliton systems. Eqs. (15) and 
( |l9| ) correspond to one of the evolution equations dis- 
cussed by Ablowitz and Ladik M. The soliton systems 



related to Eq.(19) were widely discussed in the self dual 
network equations of nonlinear inductors and capacitors 
by Wadati , Hirota and Satsuma . Wadati showed 
the corresponding soliton system to Eq. jl9| ) was derived 
from Lotka-Volterra system by Backlund transformation 
p0| . These soliton systems reduce to MKdV equation. 
Our system of difference equation with delay may be re- 
lated to soliton systems. 
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